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Beschreibung 



Rechnergestiitztes Modellierverf ahren fur das Verhalten eines 
St ahlvo lumens itiit einer VolxiiuenoberflSche 

5 

Die vorliegende Erfindung betrifft ein rechnergestiitztes Mo- 
dellierverf ahren fiir das Verhalten eines Stahlvolumens mit 
einer Volumenoberf lacher 

- wobei ein Rechner anhand eines momentanen Anf angszustands 
10 des Stahlvolumens und mindestens einer iiber die Volumen- 
oberf lache auf das St ahlvo lumen einwirkenden momentanen 
Einf lussgroBe durch Losen einer Warmeleitungsgleichung und 
einer Phasenumwandlungsgleichung einen Folgezustand des 
Stahlvolumens ermittelt , 

15 - wobei die mindestens eine Einf lussgroBe fiir eine Anzahl von 
Flachenelementen der Volumenoberf lache jeweils einen loka- 
len Einfluss umfasst und die lokalen Einfliisse iiber das je- 
weilige Flachenelement auf das Stahlvolximen einwirken, 

- wobei der Anf angszustand und der Folgezustand fiir eine An- 
20 zahl von Volumenelementen des Stahlvolumens jeweils lokale 

Anteile von modellierten Phasen des Stahls und eine einen 
lokalen Energieinhalt des Stahls beschreibende GroBe xamfas- 
sen, 

- wobei die modellierten Phasen des Stahls Austenit und eine 
25 erste weitere Phase umfassen, in die Austenit umwandelbar 

ist und die in Austenit umwandelbar ist* 

Ein derartiges Modellierverf ahren ist z. aus der DE-A-101 
29 565 bekannt, Bei dieser Schrift wird inSbesondere erstmals 

30 versucht/ die Fouriersche Warmeleitiongsgleichung selbst und 
nicht eine unkorrekte Abwandlung dieser Warmeleitungsglei- 
chung zu losen, um das thermodynamische Verhalten eines 
Stahlbandes korrekt zu beschreiben, Diese Schrift wird hier- 
mit durch Verweis auf sie in den Of fenbarungsgehalt der vor- 

35 liegenden Erfindung mit aufgenommen. 
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In der alteren, zum Anmeldezeitpunkt der vorliegenden Erfin- 
dung nicht vorver5f fentlichten deutschen Patentanmeldung 102 
51 716.9 ist ebenfalls ein derartiges Modellierverfahren be- 
schrieben, Bei diesem Modellierverfahren wird versucht, die 
5 Phasenumwandlung des Stahls anhand der Gibbsschen freien 

Enthalpien des Stahls zu modellieren . Auch diese Schrift wird 
hiermit durch Verweis auf sie in den Of f enbarungsgehalt der 

■ 

vorliegenden Erfindung mit auf genommen . 

.1 

10 Weiterhin ist ein derartiges Modellierverfahren z. B. aus dem 
Aufsatz ^^Numerische Simulation der Warmeleitung in Stahlble- 
Chen - Mathematik hilft bei der Steuerung von Kiihlstrecken" 
von W, Borchers et al,^ verof fentlicht im Unikurier der 
■ Friedrich-Alexander-Universitat Erlangen-Nurnberg, Heft 102, 

15 Oktober 2001, 27. Jahrgang, bekannt. 

Schliefilich sind noch die traditionellen Ansatze gemaB der 
Scheilschen Kegel, nach Johns on-Mehl-Avr ami und Brimacombe 
bekannt . 

20 

Bereits die exakte Modellierung des zeitlichen Terr^eraturver- 
laufs von Stahl bei der Abkiihlung, insbesondere von Stahlban- 
dern, ist entscheidend fiir die Steuerung der erf orderlichen 
Wasser- bzw. Kiihlmittelmengen einer Kiihlstrecke fiir Stahl. 

25 Denn die beim Abkiihlen auftretende Stahlumwandlung beein- 
flusst entscheidend das thermische Verhalten des Stahls bei 
der Abkiihlung, Auch werden durch den Abkiihlprozess wesentli- 
che Materialeigenschaften des Stahls beeinflusst. Da die Ab- 
kiihlung aber nicht im thermodynamischen Gleichgewicht er- 

30 folgt, ist es nicht moglich, die Umwandlung einfach durch ge- 
eignete Anpassung der WarmekapazitSt zu beschreiben. Daher 
ist zusatzlich eine exakte Modellierung der Phasenumwandlung 
des Stahls erf orderlich, urn die Kiihlstrecke korrekt steuern 
zu konnen. 

35 

Die traditionellen Ansatze des Standes der Technik arbeiten 
in der Praxis nicht in alien Fallen fehlerfrei. Insbesondere 
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weisen sie eine Reihe von systematischen Nachteilen auf . Zxim 
Ersten muss jedes Material separat parametriert werdeh, In- 
terpolationen zwischen verschiedenen Materialien sind nicht 
Oder zumindest nur eingeschrankt moglich. Ziam Zweiten werden 
5 bei den traditionellen Verfahren des Standes der Technik nur 
zwei Phasen betrachtet. Eine Erweiterung auf mehr als zwei 
Phasen ist systembedingt nicht moglich. Zum Dritten liefern 
die Verfahren des traditionellen Standes der Technik nur bei 
einer vollstSndigen Umwandlung des betrachteten Metalls eine 

10 gute Ubereinstiininung zwischen Modell und Realitat. Zum Vier- 
ten liefern die Verfahren des traditionellen Standes der 
Technik keine Information uber die bei der Phaseniomwandlung 
freigesetzte Umwandlungswarme • Die Kenntnis der Umwandlungs- 
warme ist fiir ein korrektes Losen der Warmeleitungsgleichung 

15 aber unabdingbar. 

Die Vorgehensweisen gemaS der DE-A-101 29 565 und dem Fach- 
aufsatz ^^Numerische Simulation stellen. dem gegenuber be- 

reits einen erheblichen Fortschritt dar, da sie zxunindest die 
20 Warmeleitung vollig korrekt beschreiben. Die altere deutsche 
Patentanmeldung verbessert daruber hinaus die Modellierung 
der Phasenumwandlung . Insbesondere liefert sie bereits die 
bei der Phasenumwandlung auftretende Umwandlungswarme., Auch 
diese Verfahren sind aber noch verbesserungsfahig. 

25 

Die Aufgabe der vorliegenden Erfindung besteht darin, ein Mo- 
dellierverf ahren ffir ein Metall zu schaffen, das bessere Mo- 
dellergebnisse liefert, 

30 Die Aufgabe wird dadurch gelost, 

- dass der Anf angszustand und der Folgezustand fiir mindestens 
eines der Volumenelemente auch eine lokale Konzentrations- 
verteilung mindestens eines im Stahl beweglichen Legie- 
rungselements umfassen/ 

35 - dass im Rahmen der Dmwandlungsgleichung fiir das mindestens 
eine Volvimenelement ermittelt wird, welche Konzentrationen 
des mindestens einen beweglichen Legierungselements beid- 
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seits einer ersten Phasengrenze zwischen Austenit und der 
ersten weiteren Phase vorliegen, 

- dass durch Losen eines ersten Stefan-Problems ermittelt 
wird/ ob und wie sich die Konzentrationsverteilung des min- 

5 destens einen beweglichen Legierungselements im austeniti- 
schen Bereich des betrachteten Volumenelements andert und 
Ob und um welches Ausmafi sich die erste Phasengrenze da- 
durch verschiebt, und 

- dass die lokalen Anteile der Phasen anhand einer durch das 
10 AusmaB der Verschiebung der ersten Phasengrenze bestimmten 

Lage der ersten Phasengrenze ermittelt werden* 

Das in Stahl bewegliche Legierungselement ist in der Regel 
Kohlenstoff . Alternativ oder zusatzlich ist es aber auch mog- 
15 lich, dass das Legierungselement Stickstoff ist. Es kann auch 
ein anderes Legierungselement sein, das in Stahl bevorzugt an 
Zwischengitterplatzen angeordnet ist, 

Bereits dieses Modellierverf ahren ermoglicht gegenuber dem 
20 bisher bekannten Stand der Technik einen deutlichen Fort- 
schritt. Denn mit dieser Vorgehensweise kann. - je nach dem, 
ob die weitere Phase Ferrit oder Zementit ist - das Umwand- 
lungs verhalt en von Austenit in Ferrit bzw. von Austenit in 
Zementit und ximgekehrt sehr realitatsgetreu modelliert wer- 
25 den, 

Vorzugsweise umf assen die modellierten Phasen des Stahls auch 
eine zweite weitere Phase, in die Austenit umwandelbar ist 
und die in Austenit umwandelbar ist. Denn dann ist es insbe- 
30 sondere moglich, 

- dass fiir das betrachtete Volumenelement im Rahmen der Um- 
wandlungsgleichung auch ermittelt wird, welche Konzentrati- 
onen des mindestens einen beweglichen Legierungselements 
beidseits einer zweiten Phasengrenze zwischen Austenit und 

35 der zweiten weiteren Phase vorliegen, 

- dass durch zusatzliches Losen eines zweiten Stefan-Problems 
ermittelt wird, ob und wie sich die Konzentrationsvertei- 
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lung des mindestens einen beweglichen Legierungselements im 
austenitischen Bereich des betrachteten Volumenelements an- 
dert xand ob und urn welches AusmaB sich die zweite Phasen- 
grenze dadurch verschiebt, 
5 - dass die Stef an-Probleme miteinander gekoppelt sind/ 

- dass den Phasengrenzen Flachenmafie zugeordnet sind, 

- dass ein Anteil des der zweiten Phasengrenze zugeordneten 
FlachenmaBes an der Summe der Flachenmafie bestimmt wird und 

- dass die lokalen Anteile der Phasen auch vom Anteil des der 
10 zweiten Phasengrenze zugeordneten Flachenmafies an der Summe 

der Flachenmafie abhangen. 

Mit dieser Vorgehensweise ist nicht nur die Umwandlung zwi- 
schen Austenit auf der einen Seite und Ferrit oder Zementit 

15 auf der anderen Seite sehr realitatsgetreu beschreibbar, son- 
dern insbesondere auch die Umwandlung von Austenit in Perlit 
und umgekehrt. Es ist somit das gesamte wesentliche Umwand- 
lungsverhalten von Stahl in festem Zustand korrekt beschreib- 
bar. Ferner ist mit dieser Vorgehensweise ermittelbar, ob 

20 Perlit gebildet wird oder nicht. 

Der Anteil des der zweiten Phasengrenze zugeordneten Flachen- 
maJ3es an der Summe der FlachenmaJie kann derart bestimmt wer- 
den, dass die Phasengrenzen stets nebeneinander angeordnet 
25 bledLben. Alternativ ist es aber auch moglich, diesen Anteil 
derart nachzufiihren, dass die Phasengrenzen aufeinander zu 
streben, Anhand dieses Anteils kann dann auch abgeleitet wer- 
den, ob Austenit nur in die erste weitere Phase, nur in die 
zweite weitere Phase oder sowohl in die erste als auch in die 

* 

30 zweite weitere Phase umgewandelt wird- 

Prinzipiell ist es moglich, die Ausbildung der ersten Phasen- 
grenze dreidimensional zu ermitteln. Dies fordert aber einen 
erheblichen Rechenaufwand. Vorzugsweise wird daher die Rech- 
35 nung eindimensional vorgenommen. Dies ist einfacher moglich, 
wenn das Volumenelement ein raumfiillender aristotelescher 
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Korper/ insbesondere ein Quader/ ist. Es ist also vorzugs- 
welse so, 

- dass das betrachtete Volumenelement quaderf orinig ausgebil- 
det ist und drei Grundabmessungen aufweist, 

5 - dass die erste Phasengrenze als Rechteck mit einer ersten 
Langsseite und einer ersten Querseite ausgebildet ist und 

- dass die erste Langsseite mit einer ersten der Grundabmes- 
sungen korrespondiert, die erste Querseite parallel zu ei- 
ner zweiten der Grundabmessungen verlauft und Verschiebun- 

10 gen der ersten Phasengrenze parallel zu der dritten der 
Grundabmessungen erfolgen. 

Bei Berucksichtigung auch der zweiten weiteren Phase ist es 
dann so^ 

15 - dass die zweite Phasengrenze als Rechteck mit einer zweiten 
Langsseite und einer zweiten Querseite ausgebildet ist und 

- dass die zweite Langsseite mit der ersten der Grundabmes- 
sungen korrespondiertf die zweite Querseite parallel zur 
zweiten der Grundabmessungen verlSuft und Verschiebungen 

20 der zweiten Phasengrenze parallel zur dritten der Grundab- 
messungen erfolgen* 

Wenn die Summe der Querseiten der Phasengrenzen ungefahr 
gleich dem 1,5- bis 3-fachen eines kritischen Lamellenabstan- 

25 des ist/ bei dem eine Energiebilanz, die einer seits mit dem 
Verschieben der Phasengrenzen korrespondierende Phasenumwand- 
lungen des Stahls und andererseits mit dem Verschieben der 
Phasengrenzen korrespondierende Anderungen der Flache einer 
Grenzschicht zwischen der ersten und der zweiten weiteren 

30 Phase beriicksichtigt, ausgeglichen ist, kann mittels des Mo- 
dells auch der Lamellenabstand des Perlits ermittelt werden. 
Der Lamellenabstand (das heilJt die Summe der Querseiten der 
Phasengrenzen) kann dabei insbesondere in etwa gleich dem 
Doppelten des kritischen Lamellenabsteunds sein. 

35 

Trotz der nur eindimensionalen Rechnung liefert das erfin- 
dungsgemaBe Modellierverfahren sehr realitatsgetreue Ergeb- 
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nisse, wenn der Anteil an Austenit anhand einer nicht linea- 
ren Funktion der Lage der Phasengrenze bzw« der Phasengrenzen 
ermittelt wird. 

5 Wenn die Ermittlung, in welchen Konzentrationen das mindes- 
tens eine bewegliche Legieriingselement beidseits der ersten 
Phasengrenze bzw. beidseits der ersten und beidseits der 
zweiten Phasengrenze vorliegt, anhand der Gibbsschen freien 
Enthalpien der Phasen erfolgt, ist die Ermittlung der Kon- 
10 zentrationen des mindestens einen beweglichen Legierungsele- 
ments an den Phasengrenzen besonders einfach, 

Es ist moglich, stets die gekoppelten Stef an-Probleme zu 15- 
sen und daraus abzuleiten, welche Phasen in welchem Umfang 

15 gebildet werden. In manchen Fallen ist dies sogar unumgSng- 
lich. Manchmal ist es aber mdglich, vor dem Losen der Phasen- 
umwandlungsgleichung anhand der im Anf angszustand bereits 
vorhandenen Phasen und anhand der Gibbsschen freien Enthal- 
pien vorab zu ermitteln, ob sowohl Austenit als auch die 

20 erste weitere Phase vorliegen bzw. ob zusatzlich zu Austenit 
und der ersten we,iteren Phase auch die zweite weitere Phase 
vorliegt . 

Im Einzelfall kann die Anzahl von Volumenelementen klein 
25 sein. Im Extremfall ist sie gleich Eins. In aller Regel um- 
f asst das Stahlvolumen aber eine Vielzahl von Volumenelemen- 
ten. Wenn das Stef an -Problem bzw. die Stef an-Probleme nur fur 
einen Teil der Volumenelemente gelost wird bzw. werden und 
die lokalen Anteile der Phasen der anderen Volumenelemente 
30 anhand der lokalen Anteile der Phasen des Teils der Volumen- 
elemente ermittelt werden, kann daher die zum Modellieren des 
Verhaltens des Stahlvolumens erf orderliche Rechenleistung 
deutlich reduziert werden, ohne die Aussagekraft der Modell- 
rechnung wesentlich zu verschlechtern. Die WSrmeleitungsglei- 
35 chung hingegen wird in aller Regel individuell fiir jedes Vo- 
lumenelement gelost. 
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Das Modellierverfahren kann alternativ online und in Echtzeit 
Oder aber offline ausgefuhrt warden. 

Es ist z, moglich, 
5 - dass dem Rechner ein Erstzustand und zumindest eine ge- 
wiinschte Endgr65e vorgegeben werden, 
~ dass das obige Modellierverfahren i1:erativ angewendet wird, 

- dass der Anfangszustand der ersten Iteration dem Erstzu- 
stand und der Anf aucigszustand jeder weiteren Iteration dem 

10 unmittelbar zuvor ermittelten Folgezustand entspricht und 

- dass anliand des nach einer letzten Iteration ermittelten 
Folgezustands eine erwartete Endgrofie ermittelt und mit der 
gevninschten EndgroBe verglicben wird. 

15 In diesem Fall kann das erf indungsgemaBe Modellierverfahren 
alternativ online und in Echtzeit oder offline ausgefuhrt 
werden . 

Insbesondere bei einer of fline-Ausfiihrung ist es dabei mog- 
20 lich, dass die Einf lussgroBen der Iterationen in ihrer Ge- 

samtheit einem Einf lussgroBenverlauf entsprechen,. der Rechner 
auf Grund des Vergleichs der eirwarteten EndgroBe mit der ge- 
wunschten EndgrSBe den EinflussgrSBenverlauf variiert und 
ausgehend vom Erstzustand das obige Modellierverfahren erneut 
25 ausfiihrt, bis zumindest die erwartete EndgroBe mit der ge- 
wunschten EndgroBe korrespondiert , 

Es ist aber auch moglich, dass der Rechner anhand. einer aus 
dem Anfangszustand ermittelten AnfangsgroBe und einer ge- 

30 wunschten FolgegroBe die EinflussgroBe ermittelt und eine Be- 
einflussungseinrichtung derart ansteuert/ dass das Stahlvolu- 
men entsprechend der ermittelten EinflussgroBe beeinflusst 
wird. In diesem Fall muss das Modellierverfahren online und 
in Echtzeit ausgefuhrt werden. Das Ansteuern der Beeinflus- 

35 sungseinrichtung kann dabei alternativ sof ort oder in der 
darauff olgenden Iteration erfolgen. 
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Weitere Vorteile aind Einzelheiten ergeben sich aus der nach- 
folgenden Beschreibxing eines Aus fiihrungsbei spiels in Verbin- 
dung mit den Zeichnungen, Dabei zeigen in Prinzipdarstellung 

5 FIG 1 eine Kuhlstrecke fiir ein Stahlband, 
FIG 2 das Stahlband von FIG 1 im Detail, 

FIG 3 ein Voliimenelement in perspektivischer Darstellung, 
FIG 4 ein Ablaufdiagramm/ 
FIG 5 ein weiteres Ablaufdiagramm, 
10 FIG 6 eine Warmeleitungsgleichung und eine Phasen\jmwand- 

lungsgleichung/ 

FIG 7 eine weitere warmeleitungsgleichung und eine Phasen- 

umwandlungs glei chung , 
FIG 8 ein weiteres Ablaufdiagramm, 
15 FIG 9 die Gibbsschen freien Enthalpien der Phasen des 

Stahls als Funktion des Kohlenstof f anteils bei einer 

erst en Temper at ur, 
FIG 10 die Gibbsschen freien Enthalpien als Funktion des 

Kohlenstoffanteils des Stahls bei einer zweiten Tem- 
20 peratur, 

FIG 11 ein Volumenelement in perspektivischer Darstellung, 
FIG 12 das Volumenelement von FIG 11 in der Seitenansicht, 
FIG 13 ein Volumenelement analog dem von FIG 12, 
FIG 14 ein Volumenelement in perspektivischer Darstellung, 
25 FIG 15 das Volumenelement von FIG 14 in der Draufsicht, 

FIG 16 das Volumenelement von FIG 15 im Schnitt gemaB der 

Linie XVI - XVI von FIG 15, 
FIG 17 das Volumenelement von FIG 15 im Schnitt gemafi der 

Linie XVII - XVII von FIG 15, 
30 FIG 18 ein weiteres Ablaufdiagramm und 
FIG 19 ein Detail von FIG 15. 

GemaB FIG 1 soil z. B, ein Stahlband 1 gekiahlt werden, damit 
es vorgegebene Endeigenschaften annimmt. Hierzu wird das 
35 Stahlband 1 z. B. in einer Kuhlstrecke 2 mit einem Kiihlmedium 
3 - in der Kegel Wasser - beaufschlagt . 
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Das Stahlband 1 weist eine Bandbreite b und eine Banddicke d 
auf . Bs weist ferner eine - prinzipieli beliebige - Bandlange 
auf . Es durchlauft die Kiihlstrecke 2 mit einer Bandgeschwin- 
digkeit v, 

5 

Die Bandgeschwindigkeit v wird kontinuierlich erfasst und ei- 
nem Rechner 4 zugefiibrtr der die Kuhlstrecke 2 steuert. Der 
Rechner 4 ist somit in der Lage, in an sich bekannter Weise 
eine Wegverfolgung des Stahlbandes 1 zu realisieren, Dem 
10 Rechner 4 ist also bekannt, auf welche Bereiche des Stahlban- 
des 1 eingewirkt wird, wenn das Stahlband 1 mittels Aufbrin- 
gungseinrichtungen 5 der Kuhlstrecke 2 mit dem Kiihlmediiam 3 
beaufschlagt wird. 

15 Der Rechner 4 ist mit einem Computerprogramm 6 programmiert, 
das dem Rechner 4 uber einen Datentrager 7, z. B, eine CD-ROM 
If zugefiihrt wird. Auf dem Datentrager 7 ist das Computerpro- 
gramm 6 in (ausschlieiilich) maschinenlesbarer Form gespei- 
chert. Das Computerprogramm 6 wird vom Rechner 4 iibernommen 

20 und in einem Massenspeicher 8/ z. B. einer Festplatte 8, des 
Rechners 4 hinterlegt. Dadurch ist der Rechner 4 in der Lage, 
bei Aufruf des Computerprogramms 6 ein Modellierverf ahren fur 
das Stahlband 1 bzw. fur dessen einzelne Bereiche (= das 
Stahlvolumen 1) auszufiihren, das nachfolgend naher beschrie- 

25 ben wird. 

GemSfi FIG 2 wird hierzu das Stahlvolxomen 1 innerhalb des 
Rechners 4 in Volumenelemente 9 zerlegt. Soweit ein Volumen- 
element 9 dabei nicht allseitig von anderen Volumenelementen 
30 9 umgeben ist, sind dem jeweiligen Volumenelement 9 ein oder 
zwei FlSchenelemente 10 zugeordnet. Die Flachenelemente 10 
bilden in ihrer Gesamtheit eine Volumenoberf lache des Stahl- 
bandes 1 bzw. des betrachteten Bereichs des Stahlbandes 1. 

35 Der Rechner 4 realisiert auf Grund der Programmierung mit dem 
Computerprogramm 6 insbesondere ein Modell, in dem fur jedes 
Volumenelement 9 die thermischen Kopplungen mit seiner Umge- 

10 
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bung beriicksichtigt sind. Jedes Voliimenelement 9 ist dabei - 
siehe FIG 3 - quaderformig ausgebildet. Es weist also drei 
Grundabmessungen A, C auf , die in der Kegel parallel zur 
Richtung der Bandgeschwindigkeit v, zur Bandbreite b und zur 
5 Banddicke d orientiert sind. 

In einer ersten Ausgestaltung des erf indungsgemaBen Model- 
lierverf ahrens wird dem Rechner 4 gemaB FIG 4 in einem 
Schritt SI ein Erstzustand Z zugefuhrt. Der Erstzustand Z urn- 
10 fasst fiir jedes Volumenelement 9 des Stahlbandes 1 zunachst 
lokaie Anteile pi, p2, p3 von modellierten Phasen des Stahls, 
Die Phasen kSnnen insbesondere Ferrit (Anteil pi) / Zementit 
(Anteil p2) und Austenit (Anteil p3) sein, 

15 Der Erstzustand Z lomfasst weiterhin fur jedes Volumenelement 
9 des Stahlbandes 1 eine GrSfie H, die einen lokalen Energie- 
inhalt des Stahls des Volumenelements 9 beschreibt, Bei- 
spielsweise kann diese GroBe H die Enthalpie H des Volumen- 
elements 9 sein. Alternativ kommen auch die Temperatur oder 

20 die Entropie in Frage. 

SchlieBlich umfasst der Erstzustand Z fiir mindestens eines 
der Volumenelement e 9 des Stahlbandes 1, vorzugsweise fur je- 
des Volumenelement 9 des Stahlbandes 1, noch eine lokaie Kon- 
25 zentrationsverteilxing K mindestens eines im Stahl beweglichen 
Legierungselements • Das bewegliche Legierungselement kann 
insbesondere Kohlenstoff sein. Alternativ oder zusStzlich 
kommt z. B. auch Stickstoff in Frage. 

30 Dem Rechner 4 wird im Rahmen des Schrittes SI weiterhin zu- 
mindest eine gewiinschte EndgroBe f* zugefuhrt. Gegebenen- 
f alls kSnnen dem Rechner 4 im Rahmen des Schrittes SI auch 
noch gewiinschte ZwischengroBen zugefiihrt werden, so dass dem 
Rechner 4 gegebenenf alls sogar ein gewiinschter zeitlicher 

35 Verlauf der Gr5Be vorgegeben werden kann. 
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GemaB FIG 4 wird nun in einem Schritt S2 zunachst ein Ein- 
.flussgroBenverlauf festgelegt. Der Einf lussgroBenverlauf inn- 
fasst fiir eine Vielzahl von zeitlich unmittelbar aufeinander- 
folgenden Zeitpunkten jeweils eine EinflussgroBe Die Ein- 
5 flussgroBe W entspricht beispielsweise der auf das Stahlvolu- 
men 1 (also das Stahlband 1 bzw. dessen betrachteten Bereich) 
auf zijODringenden Kiihlitiittelmenge unter Beriicksichtigung ande- 
rer Einflusse auf das Stahlvolumen 1/ wie beispielsweise von 
Transportrollen, der Warmekonvektion, der Warmeabstrahlung 
10 usw.. Sie umfasst fiir die Vielzahl von Flachenelementen 10 
des St ahlvo lumens 1 (siehe FIG 2) jeweils einen lokalen Ein- 
fluss. Die lokalen Einflusse wirken dann iiber das jeweilige 
Flachenelement 10 "auf das Stahlvolumen 1 ein. 

15 Gemafi FIG 4 wird nunmehr in einem Schritt S3 ein Anfangszu- 
stand ZA des Stahlvolumens 1 gleich dem Erstzustand Z und 
eine Zeitbasis t auf einen Anf angszeitpunkt tO gesetzt, Der 
Anfangszeitpunkt tO entspricht dabei in der Kegel dem Zeit- 
punkt, zu dem das betrachtete Stahlvolumen 1^ also z. B. ein 

20 Abschnitt des Stahlbandes 1, in die Kuhlstrecke 2 eintritt. 

Sodann wird in einem Schritt 34 anhand des Einf lussgroBenver- 
laufs die zu dem durch die Zeitbasis t bestimmten Zeitpunkt 
auf zubringende EinflussgroBe W ermittelt. In einem Schritt S5 

25 ermittelt der Rechner 4 dann anhand des momentanen Anf angszu- 
stands ZA des Stahlvolumens 1 und der iiber die Volumenober- 
flSche auf das Stahlvolumen 1 einwirkenden momentanen Ein-~ 
flussgroBe W einen Folgezustand ZF des Stahlvolumens 1. Der 
JRechner 4 lost dabei eine Warmeleitungsgleichung und eine 

30 Phasenumwandlungsgleichung. Der Folgezustand ZF umfasst die- 
selben Elemente pi, p2, p3^ H wie der Anfangszustand ZA. 

Sodann wird in einem Schritt S6 die Zeitbasis t um eine zeit- 
liche Schrittweite At inkrementiert . Als nachstes wird in 
35 einem Schritt S7 gepruft, ob die Zeitbasis t einen Endzeit- 
punkt tl erreicht hat. Der Endzeitpunkt tl entspricht dabei 
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in der Kegel dem Zeitpunkt, zu dem das betrachtete Stahlvolu- 
men 1 die Kiihlstrecke 2 wieder verlasst. 

Wenn der Endzeitpunkt tl noch nicht erreicht ist, wird in ei- 
5 nem Schritt S8 der Anf angszustand ZA gleich dem unmittelbar 
zuvor ermittelten Folgezustand ZF gesetzt und sodann zxm 
Schritt S4 zuruckgespmngen . 

Wenn hingegen der Endzeitpunkt tl erreicht wurde, wird die 
10 aus den Schritten S4 bis S8 bestehende Schleife verlassen. 
Wenn danach keine weiteren Mafinahmen ergriffen werden, wird 
nur ein sogenannter Prozessbeobachter realisiert. Vorzugs- 
weise aber wird, wie in FIG 4 dargestellt, in einem Schritt 
S9 anhand des nunmehr ermittelten Folgezustands ZA eine er- 
15 wartete EndgroBe f'f z. B. die Temperatur oder die Material- 
harte, ermittelt und mit der gewunschten Endgrofie f* vergli- 
chen. Wenn die erwartete Endgrofie f mit der gewunschten End- 
grofie nicht korrespondiert, wird - die entsprechende Re- 
chenleistung vorausgesetzt - in einem Schritt 310 der Verlauf 
20 der Einflussgrofie W variiert. Sodann wird ziom Schritt S3 zu- 
ruckgesprungen . Anderenfalls ist der erf orderliche Kiihlmit- 
telmengenverlauf ermittelt, so dass der Rechner 4 nunmehr die 
Kfihlstrecke 2 in einem Schritt Sll entsprechend ansteuern 
kann- 

25 

Der Vollstandigkeit halber sei dabei erwahnt, dass der Rech- 
ner 4, wenn ihm aixBer der gewunschten Endgrofie f * auch ge- 
wiinschte ZwischengroBen vorgegeben werden, selbstverstandlich 
den Kiihlmittelmengenverlauf erst bis zu der ersten gewiinsch- 
30 ten Zwischengrofie ermittelt, dann bis zur zweiten gewunschten 
Zwischengrofie usw., bis der gesamte Kuhlmittelmengenverlauf 
bis zur gewunschten Endgrofie f* ermittelt ist. 

Gemafi den FIG 1 und 4 wird das obenstehend beschriebene Mo- 
35 dellierverf ahren online und in Echtzeit ausgefiihrt, Es ist 
aber selbstverstandlich ohne Weiteres auch eine Offline-Aus- 
fiihrung moglich. Dies ist in FIG 1 dadurch angedeutet, dass 
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die Verbindxmg der Kiihlstrecke 2 mit dem Rechner 4 als offen- 
bare Verbindung angedeutet ist, also nicht unmittelbar eine 
Ansteuerung der Kiihlstrecke 2 erfolgen muss, Weiterhin ist 
aus diesem Grund der Schritt Sll in FIG 4 nur gestrichelt 
5 eingezeichnet • 

Die Vorgehensweise gemafi FIG 4 erfordert einen sehr hohen Re- 
chenaufwand. Wenn die zur Ausfiihrung des Verfahrens gem^ FIG 
4 erforderliche Rechenlei stung nicht zur Verfiigung steht^ 
10 aber dennoch eine Online-Steuerung unter Nutzung des erfin- 
dungsgemaBen Modellierverf ahrens erfolgen soli, wird online 
und in Echtzeit ein Verfahren ausgefiihrt/ das nachfolgend in 
Verbindung mit FIG 5 naher beschrieben wird. 

15 GemaJi FIG 5 nimmt der Rechner 4 in einem Schritt S12 analog 
zum Schritt SI den Erstzustand Z entgegen. In einem Schritt 
S13 setzt der Rechner 4 analog zum Schritt S3 den Anfangszu- 
stand ZA gleich dem Erstzustand Z und die Zeitbasis t auf den 
Zeitpunkt tO. Sodann ermittelt der Rechner 4 in einem Schritt 

20 S14 eine gewunschte Folgegrofie f* Oder nimmt sie entgegen. 

In einem Schritt S15 ermittelt der Rechner 4 analog zum 
Schritt S9 anhand des Anf angszustands ZA eine AnfangsgrSfie f . 
Anhand der im Schritt S15 ermittelten AnfangsgroBe f und der 

25 gewiinschten FolgegroBe f* ermittelt der Rechner 4 dann in ei- 
nem Schritt S16 die Einflussgrofie W. In einem Schritt S17, 
der mit dem Schritt Sll korrespondiert, steuert der Rechner 4 

/ schlieBlich die Kiihlstrecke 2 entsprechend der ermittelten 

Einf lussgroBe W an. Das Ansteuern erfolgt dabei, wie in FIG 5 

30 dargestellt, vorzugsweise sofort, Gegebenenf alls kann aber 
auch ein Ansteuern im nachsten Zyklus erfolgen. 

Die nachf olgenden Schritte S18 bis S21 korrespondieren mit 
den Schritten S5 bis S8 von FIG 4. Beziiglich dieser Schritte 
35 SI 8 bis S21 wird daher von einer detaillierten Erlauterxong 
abgesehen . 
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Das Verfahren gemafi FIG 5 wird vorzugsweise z* B, bei einer 
Regelung eines die Kiihlstrecke 2 durchlaufenden Stahlbandes 1 
auf einen Verlauf der gewiinschten FolgegroBe f* eingesetzt. 
In diesem Fall muss das obenstehend in Verbindung mit FIG 5 
5 beschriebene Verfahren selbstverstandlich individuell fiir je- 
den einzelnen i^Dschnitt des Stahlbandes 1 ausgefiihrt werden. 
Auch muss fiir jeden Abschnitt jeweils diejenige der Aufbrin- 
gungseinrichtungen 5 angesteuert werden, in deren Einflussbe- 
reich sich der betrachtete Abschnitt gerade befindet. Dies 
10 wird durch die eingangs erwahnte Wegverfolgung gewahrleistet . 
Das Verfahren gemafi FIG 5 wird in diesem Fall ferner parallel 
fur alle in der Kiihlstrecke 2 befindlichen Abschnitte des 
Stahlbandes 1 ausgefiihrt. 

15 Die jeweils gewiinschte Folgegrofie f * kann dem Rechner 4 ex- 
plizit vorgegeben werden. Es ist aber auch moglich, dass dem 
Rechner 4 lediglich die gewiinschte FolgegroBe f* z. B. beim 
Austritt des Stahlbandes 1 aus der Kiihlstrecke 2 Oder beim 
Erreichen eines Haspels, auf den das Stahlband 1 aufgehaspelt 

20 wird, vorgegeben wird. In diesem Fall ermittelt der Rechner 4 
selbsttatig anband einer vorbestimmten Ermittlungsvorschrift 
die gewiinschten FolgegroBen f* fiir die einzelnen. Iterationen. 

Die Regelung gemaB FIG 5 kann selbstverstandlich auch ander- 
25 weitig strukturiert werden, z. B. so, wie sie in der alteren 
DE 103 21 792.4 beschrieben ist. 

Wie erwahnt, werden in den Schritten S5 und S18 eine Warme- 
leitungsgleichung und eine Phaseniimwandlungsgleichung gelost. 

30 

Ein Beispiel einer WSrmeleitungsgleichung und einer Phasenum- 
wandlungsgleichung sind in FIG 6 dargestellt. Dort ist ein 
dreidimensionaler Ansatz fiir die Warmeleitung getroffen. Bei 
Bandern kann dieser Ansatz in aller Regel aber auf den eindi- 
35 mensionalen Ansatz gemafi FIG 7 vereinfacht werden, da der 
Warmefluss in Bandlangsrichtung und Bandquerrichtung in der 
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Kegel vernachlassigbar ist. 3/3x ist datoei die Ortsableitiing 
in Banddickenrichtung . 

Ersichtlich sind die Warmeleitungsgleichung und die Phasenum- 
5 wandlungsgleichung miteinander gekoppelt. Zur Losiang der War- 
meleitungsgleichung und der Phaseniomwandlungsgleichung (« 
Schritt S5r S18) wird daher gemafi FIG 8 wie folgt vorgegan- 
gen: 

10 In einem Schritt S22 wird zunachst die lokale Temperatur T 

ermittelt. Die Ermittlung erfolgt anhand der Enthalpie H und 
der Anteile pi, p2, p3 der Phasen Ferrit, Zementit und Auste- 
nit. Auf Grund der geringen GroBe des betrachteten Volumen- 
elements 9 kann namlich angenommen werden, dass die lokale 

15 Temperatur T innerhalb des Volumenelements 9 konstant ist. 
Damit weisen auch die Phasen des Stahls innerhalb dieses Vo- 
lumenelements 9 dieselbe lokale Temperatur T auf. Die Enthal- 
pie H des Volxamenelements 9 lasst sich somit als 

20 = piHiiT) + piHig:) + pm3{T) ( i > 

schreiben. Da ferner fiir jede der Phasen die Enthalpien Hi, 
H2 bzw, H3 der jeweiligen Phase als Funktion der lokalen Tem- 
peratur T - gegebenenf alls bei Austenit unter Berucksichti- 
25 gung des Kohlenstof f anteils - eineindeutig bestimmt sind, 

lasst sich durch die Gleichung 1 die lokale Temperatur T ohne 
Weiteres ermitteln. 

Fur die im Schritt S22 ermittelte lokale Temperatur T werden 
30 dann in einem Schritt S23 - getrennt fiir die Phasen Ferrit, 
Zementit und Austenit - die Gibbsschen freien Enthalpien Gl 
bis G3 der einzelnen Phasen als Funktion des Anteils des be- 
weglichen Legierungselements ermittelt. Beispiele derartiger 
Verlaufe sind in den FIG 9 und 10 dargestellt. Bei dem Bei- 
35 spiel gemSB FIG 9 liegt die lokale Temperatur T dabei ober- 
halb der Umwandlungs temperatur von Austenit in Per lit, bei 
dem Beispiel gemafi FIG 10 darunter. 
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Es ist moglichr stets anzunehmen, dass alle drei Phasen Fer- 
rit, Zementit und Austenit vorliegen und erst spater - siehe 
die nachfolgenden Erlauterungen zu FIG 18 - zu entscheiden, 
welche Phasen vorliegen und welche Phasenumwandlungen erfol— 
5 gen. Gemafi FIG 8 wird aber zunachst in einem Schritt S24 die 
Anzahl der zu betrachtenden Phasen ermittelt. Dies geschieht 
wie folgt: 

Zunachst wird das System Ferrit- Austenit heraus gegriffen. 

10 Fur dieses System wird iiberpriift, ob es bei dem im Voliomen- 
element 9 vorhandenen Gesamtanteil des beweglichen Legie- 
rungselements rein ferritischr rein austenitisch oder ge- . 
mischt stabil ist und welche Phasenverteilung gegebenenfalls 
im stabilen Zustand vorliegt, Der Gesamtanteil des bewegli- 

15 Chen Legierungselements kann anhand der Konz ent rat ions vert ei- 
lung ohne Weiteres ermittelt werden. 

Die Ermittlung der stabilen Phase (n) und gegebenenfalls der 
Phasenverteilung geschieht dadurch, dass versucht wird, ein 

20 Minimum fiir die gesamte Gibbssche freie Enthalpie G eines 

derartigen Systems zu ermitteln. Die konkrete Vorgehensweise 
ergibt sich dabei so, wie sie in den Seiten 16 bis 18 der al- 
teren deutschen Patentanmeldung 102 51 716.9 beschrieben ist. 
Diese Vorgehensweise liefert zugleich auch Konzentrationen kl 

25 und k3, in denen das bewegliche Legierungselement - typi- 
scherweise Kohlenstoff - an einer etwaigen Phasengrenze 11 
zwischen Ferrit und Austenit im ferritischen bzw. im austeni- 
tischen Bereich gegebenenfalls vorliegt. 

30 In analoger Weise werden sodann die Mischungssysteme Ferrit- 
Zementit (= Perlit) und Austenit-Zementit untersucht. Die Un- 
tersuchung des Mischungssystems Austenit-Zementit liefert da- 
bei zugleich auch Konzentrationen k2 und k4, in denen das be- 
wegliche Legierungselement an einer etwaigen Phasengrenze 12 

35 zwischen Zementit und Austenit im zementitischen bzw. im 
austenitischen Bereich gegebenenfalls vorliegt. 
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Die drei fiir die drei Zweiphasen-Systeme ermittelten miniina- 
len Gibbsschen freien Enthalpien G warden miteinander vergli- 
chen und das Zweiphasen-System mit der insgesamt minimalen 
Gibbsschen freien Enthalpie G herausgegrif fen. Wenn dieses 
5 Zweiphasen-System Austenit enthSlt, lasst sich anhand der 
Gibbsschen freien Enthalpien ermitteln, ob der stabile Zu- 
stand eine oder zwei Phasen umfasst, welche Phase dies gege- 
benenfalls ist bzw, welche Phasen dies gegebenenf alls sind 
und welche Konzentrationen kl, k3 bzw. k2^ k4 des beweglichen 
10 Legierungselements beidseits der Phasengrenze 11 von Austenit 
zu Ferrit bzw. der Phasengrenze 12 von Austenit zu Zementit 
jeweils vorliegen. 

Aus dem Anf angszustand ZA ist ferner bekannt, welche Phasen 
15 anfanglich zu welchen Anteilen pi bis p3 vorliegen. 

Wenn das Zweiphasen-System mit der minimalen Gibbsschen 
freien Enthalpie G Austenit enthalt, lasst sich also durch 
eine Zusainmenschau der anfSnglich vorliegenden Phasen und der 

20 als stabil ermittelten Phasen ermitteln^ ob ein Reinsystem 
mit nur einer der Phasen Ferrit, Austenit und Zementit oder 
ein Mischsystem mit zwei oder sogar drei der Phasen Ferrit, 
Austenit und Zementit betrachtet werden muss. Auf Grund des 
UmstandS/ dass anhand der Gibbsschen freien Enthalpien G, Gl 

25 bis G3 auch ermittelbar ist, welche Phase stabil ist bzw. 

welche Phasen stabil sind, ist dartiber hinaus auch die Rich- 
tung einer etwaigen Phasenumwandlung bereits bekannt. 

Wenn das Zweiphasen-System mit der insgesamt minimalen 
30 Gibbsschen freien Enthalpie G das Perlitsystem ist und der 
Anf angszustand ZA kein Austenit enthalt (p3 = 0) , liegt ein 
vollstandig umgewandeltes Gefiige, also ein Zweiphasen-System 
Ferrit- Zementit vor, das keiner Phasenxamwandlung mehr unter- 
liegt. Wenn hingegen der Anf angszustand ZA Austenit enthalt 
35 (p3 > 0) und das Perlitsystem die mnimale Gibbssche freie 
Enthalpie G aufweist, lasst sich nicht ohne Weiteres ent- 
scheiden, ob nun eine oder beide der Phasen Ferrit und Zemen- 
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tit gebildet werden. In diesem Fall wird daher vorlaufig an- 
genommen, dass beide Phasen gebildet werden, also ein Drei- 
phasensystem Ferrit-Zementit-Austenit zu betrachten ist- 

5 In einem Schritt S25 wird als nachstes gepriift, ob in dem be- 
trachteten Volumenelement 9 mehr als eine Phase vorliegt. 
Wenn dies nicht der Fall ist, erfolgt selbstverstandlich 
keine Phasenuntwandlung. Dennoch wird in einem Schritt S26 fur 
dieses Voliamenelement 9 eine Dif fusionsgleichung fiir das be- 

10 wegliche Legierungs element angesetzt und gelost, "um etwaige 

Konzentrationsschwankungen des beweglichen Legierungselements 
innerhalb des betrachteten Volumenelements 9 auszugleichen- 
Diese Vorgehensweise ist insbesondere bei Austenit, bei dem 
der Gehalt an Kohlenstoff stark schwanken kann, von Bedeu- 

15 tung. Der Ansatz und die Losung einer derartigen Dif fusions- 
gleichung ist Fachleuten allgemein bekannt. Auf den Ansatz 
der Dif fusionsgleichung und deren Losung wird daher nachste- 
hend nicht naher eingegangen, 

20 Wenn in dem betrachteten Volumenelement 9 mehr als eine Phase 
vorliegt, wird als nachstes in einem Schritt S27 abgepriift, 
ob alle drei Phasen vorliegen. Wenn dies nicht der Fall ist, 
liegt ein Zweiphasen-System vor. In diesem Fall wird in einem 
Schritt 528 iiberpriift, ob eine der beiden vorliegenden Phasen 

25 Austenit ist. Wenn dies nicht der Fall ist, liegt ein voll- 
standiges Perlit-System vor, also eine Schichtstruktur aus 
Ferrit und Zementit. Diese Struktur ist im Wesentlichen sta- 
bil. In diesem Fall miissen daher keine weiteren MaJBnahmen er- 
griffen werden. 

30 

Wenn hingegen eine der Phasen Austenit ist und somit die an- 
dere Phase Ferrit oder Zementit, existiert eine Phasengrenze 
11 zwischen Austenit und Ferrit (siehe FIG 11 und 12) oder 
eine Phasengrenze 12 zwischen Austenit und Zementit (siehe 
35 FIG 11 und 13) • In diesem Fall wird als vereinfachende Mafi- . 
nahme angenommen, dass die Phasengrenze 11 bzw. 12 als Recht- 
eck ausgebildet ist, das eine Langsseite und eine Querseite 



19 



wo 2005/076092 PCT/EP2004/053709 



aufweist. Die Langsseite und die Querseite des Rechtecks kor- 
respondieren dabei iriit den Grundabmessungen A und B des be- 
trachteten Volumenelements 9. Verschiebungen der Phasengrenze 
11, 12 erfolgen in diesem Fall, wie insbesondere aus FIG 11 
5 ersichtlich ist, parallel zur Grundabmessung C. 

Wie Fachleuten allgemein bekannt ist und auch aus FIG 12 er- 
sichtlich ist, ist in Ferrit nur sehr wenig Kohlenstoff ent- 
halten. An der Phasengrenze 11 zum austenitischen Bereich 
10 hingegen erfolgt ein Konzentrationssprung. Die Konzent ratio- 
nen kl, k3, in denen das beweglichie Legierungselement an der 
Phasengrenze 11 vorliegt, sind dabei die Konzentrationen kl, 
k3, die obenstehend im Schritt S24 fiir das Ferrit-Austenit- 
System ermittelt wurden. 

15 

Wenn, wie in FIG 12 dargestellt, eine Umwandlung von Austenit 
in Ferrit erfolgen soil, muss der ^^Konzentrationsberg'' des 
beweglichen Legierungselement s im austenitischen Bereich in 
der NShe der Phasengrenze 11 abdif fundieren. Umgekehrt muss, 

20 wenn eine Umwandlung von Ferrit in Austenit erfolgen soli, 

das ''^Konzentrationstal" des beweglichen Legierimgselements im 
austenitischen Bereich in der NShe der Phasengrenze 11 stan- 
dig aufgefullt werden. Die numerische Oder analytische Losung 
dieser Aufgabe ist allgemein als Stefan-Problem bekannt. Sie 

25 erfolgt in einem Schritt S29. Das Stefan -Problem lasst sich 
im vorliegenden Fall wie folgt formulieren: 

30 D ist dabei die - gegebenenf alls temperaturabhangige - Dif fu- 
sionskonstante des beweglichen Legierungselement s im Auste- 
nit. K' ist die Konzentration K. Die Indizes t und x bedeuten 
die Abieitung nach der Zeit bzw. dem Ort in Verschieberich- 
tung der Phasengrenze 11. Die Verschieberichtung der Phasen- 

35 grenze 11 muss dabei nicht notwendigerweise die Banddicken- 
richtung sein. 
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Fur die Gleichiing 2 sind die Anf angsbedingungen durch die 
ortliche Konzentrationsverteilung K' des beweglichen Legie- 
rungselements im Anf angszustand ZA definiert. Beim Losen des 
Stefan-Problems miissen dann nur noch die Randbedingungen be- 
5 trachtet werden, dass im ferristischen bzw* im austenitischen 
Bereich an der Phasengrenze 11 die Konzentrationen kl bzw. k3 
vorliegen, dass das bewegliche Legierungselement das betref- 
fende Volumenelement 9 nicht verlassen kann und dass sich die 
Verschiebung 5x' der Phasengrenze 11 gemaB der Stefan-Bedin- 
10 gung zu 

(k3-kl)5xWbt = -DK\ (3) 

ergibt, wobei 5t die beim Losen des Stefan-Problems benutzte 
15 zeitliche Schrittweite ist. Diese kann dabei ein Bruchteil 

{l/2r 1/3, 1/4 ..•) der zeitlichen Schrittweite At Oder aber 
gleich der zeitlichen Schrittweite At sein. 

Durch Losen des Stefan-Problems im Schritt S2 9 wird also er-. 
20 mittelt, ob und wie sich die Konzentrationsverteilung K bzw* 
K' des beweglichen Legierungselement s im austenitischen Be- 
reich des betrachteten Volumenelement s 9 andert. Zugleich 
wird dabei auch ermittelt, ob und um welches AusmaB 5x' sich 
die Phasengrenze 11 dadurch verschiebt. 

25 

Falls die weitere Phase nicht Ferrit, sondern Zementit ist, 
ergibt sich ein prinzipiell gleich gearteter LSsungsansatz 
und auch eine prinzipiell gleich geartete Losung. Der Unter- 
schied ist nur, dass bei der Bildung von Zementit im Zementit 

30 Kohlenstoff stark angereichert ist und im austenitischen Be- 
reich in der Nahe der Phasengrenze 12 zum Zementit Kohlen- 
stoff verarmt ist. Auch dies ist Fachleuten allgemein bekannt 
und in FIG 13 dargestellt. Es andern sich also die Randbedin- 
gungen. Die Verschiebung 5x" der Phasengrenze 11 gemaB der 

35 Stef an-Bedingung ergibt sich in diesem Fall zu 

(fc4- ife2)&c75r = -DK'^ ( 4 ) 



21 



wo 2005/076092 PCT/EP2004/053709 



K" ist dabei wieder die Konzentration K. 

Wenn hingegen im Schritt S28 festgestellt vmrde, dass alle 
drei Phasen vorliegen bzw. vorliegen kSnnten^ ergibt sich ein 
5 koinplizierteres Problem. Denn dann muss der Ansatz gewahlt 
werden, dass eine Schichtstruktur existiert^ die - siehe FIG 
14 und 15 - abwechselnd aus einer Ferritschicht 13 und einer 
Zementitschicht 14 besteht. Die Schichtstruktur grenzt an ei- 
nen austenitischen Bereich 15 an. Die FIG 16 und 17 zeigen 
10 jeweils einen Verlauf der Konzentration K' bzw, K" des beweg- 
lichen Legierungselements in einer Ferritschicht 13 bzw. ei- 
ner Zementitschicht 14 und im austenitischen Bereich 15 vor 
diesen Schichten 13, 14. 

* 

15 Im Rahmen des Modells wird dabei nach wie vor angenommen, 

dass die Phasengrenzen 11, 12 als Rechtecke ausgebildet sind, 
die jeweils eine Langsseite land eine Querseite aufweisen. Die 
Langsseiten korrespondieren weiterhin mit der ersten Grundab- 
messung A. Die Querseiten verlauf en nach wie vor parallel zur 

20 zweiten Grundabmessung B. Auch wird nach wie vor angenommen, 
dasjs Verschiebungen der Phasengrenzen 11, 12 parallel zur 
dritten Grundabmessung C erfolgen. 

Fiir jede der Phasengrenzen 11 und 12 wird nunmehr in einem 
25 Schritt S30 ein Stefan-Problem angesetzt und gelost. Der 
Schritt 330 ist in FIG 18 detaillierter dargestellt. 

Das Stefan-Problem fur die Phasengrenze 11 zwischen Ferrit 
und Austenit gehorcht gemaB FIG 18, siehe dort Schritt S31, 
30 folgender Gesetzmafiigkeit : 

K' ist dabei die Konzentration des beweglichen Legierungsele- 
35 ments vor der Phasengrenze 11. 
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Das Stefan-Problem fiir die Phasengrenze 12 zwischen Zementit 
und Austenit gehorcht folgender GesetzitiaJiigkeit : 



(6) 



ist dabei - analog zu - die Konzentration des bewegli- 
chen Legierungselements vor der Phasengrenze 12, 

LI und L2 Bind Kopplungsterme • Sie sind - siehe FIG 14 — 
10 Funktionen eines Lamellenabstands 1 der Schichtstruktur und 
eines Anteils q der Zementitphase an der Schichtstruktur so- 
wie dessen zeitlicher Ableitung. Ausgeschrieben kSnnen sie z 
B. lauten: 



15 



11 = 



1 



und 



(7) 



L2 = 



1 



I * 



(8) 



Wie aus den Gleichiangen 5 bis 8 ersichtlich ist, sind die 

r 

20 beiden Stef an-Probleme auf Grund der Kopplungsterme LI, L2 

auf der rechten Seite der Gleichxingen 5 und 6 miteinander ge- 
koppelt . 

Auch in den Gleichungen 5 bis 8 stehen die Indizes t und x 
25 wieder fair die Ableitung nach der Zeit bzw. dem Ort in Ver- 
schieberichtung der Phasengrenze 11 bzw. 12. 



Fur die Verschiebungen 5x' , Sx" der Phasengrenzen 11, 12 gel- 
ten weiterhin die Stef an-Bedingungen (siehe Gleichungen 3 und 
30 4) . Es ist daher moglich, die Verschiebungen 5x' , 5x" gemafi 
den Gleichungen 3 und 4 in einem Schritt S32 gleichzusetzen. 
Dann ist der Anteil q derart zia bestiiranen, dass die Verschie- 
bung 8x' der Phasengrenze 11 gemafi Gleichung 3 und die Ver- 
schiebung 5x" der Phasengrenze 12 geitiMB Gleichung 4 denselben 
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Wert annehmen. Der Anteil q ist in diesem Fall also derart zu 
bestimmen, dass die Phasengrenzen 11, 12 stets nebeneinander 



5 Die Gleichheit der Verschiebungen 5x' , 6x" der Phasengrenzen 
11, 12 kann also dadurch erreicht werden, dass der Anteil q 
der Zement it phase an der Schichtstruktur des Perlit geeignet 
gewahlt wird. Man kann also - z. B. durch Einsetzen und Aus- 
rechnen (trial and error) - ermitteln, bei welchem Anteil q 
10 der Zementitphase die Verschiebungen 5x' , Sx'' der Phasengren- 
zen 11, 12 ubereinstimmen. 

Anhand des so ermittelten Anteils q lasst sich dann auch 
bestimmen, ob tatsachlich Perlit gebildet wird Oder ob nur 
15 eine der beiden Phasen Ferrit und Zementit gebildet wird. 

Denn nur wenn der Anteil q zwischen Null und Eins liegt, wird 
Perlit gebildet. Ist der Anteil q hingegen grower als Eins, 
wird ausschliefilich Zementit gebildet. Ist er umgekehrt klei- 
ner als Null, wird ausschliefilich Ferrit gebildet. 



Der Anteil q kann bei dieser Vorgehensweise aus numerischen 
Griinden beliebige Werte annehmen, also insbesondere auch 
kleiner als Null Oder gr5£er als Eins sein. Diese Werte sind 
aber physikalisch sinnlos. Der Anteil q wird daher in Schrit- 
25 ten 833 bis 336 gegebenenf ails entsprechend korrigiert. 

Die Vorgehensweise gemafi Schritt 32 kann zu numerischen Prob- 
lemen fuhren. Alternativ ist es daher moglich, iinterschiedli- 
Che Verschiebungen 5x' , 5x" zuzulassen. In diesem Fall wird 
30 der Anteil q des Zementit an der Schichtstruktur gem^a 



angeordnet bleiben. 




(9) 



35 



in Schritten S32a und S32b derart nachgefiihrt, dass die Lagen 
der Phasengrenzen 11, 12 aufeinander zustreben. Die Schritte 
S32a und S32b werden dabei gemSB FIG 18 anstelle des Schrit- 
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tes S32 ausgefuhrt. a ist eine geeignete Proportionalitats- 
konstante, Ihr Wert ist grofier als Null. 

In diesem Fall wird eine mittlere Verschiebung 5x der Phasen- 
5 grenzen 11, 12 zu 



10 



def iniert . 



Unabhangig davon, welche dieser beiden Vorgehensweisen 
(Schritt S32 Oder Schritte S32a, S32b) ergriffen wird,. wird 
aber in beiden FS.llen erst durch die BestimmiKig des Anteils q 
ermittelt und abgeleitet, welche Phasen in welchem Urafang ge- 
15 bildet werden. 

Wenn Perlit gebildet wird, der Anteil q also zwischen Null 
und Eins liegt, muss noch der Lamellenabstand 1 bestimmt wer- 
den. Dies geschieht gemaB FIG 19 wie folgt: 

20 

Mit der Verschiebung 5x der Phasengrenzen 11, 12 und dem An- 
teil q ist bekannt, welche Phasenumwandlungen erfolgen und in 
welchen Anteilen sie erfolgen. Die Phasenumwandlungen liefern 
einen Beitrag 5e1 zur Energiebilanz . Der Beitrag 8e1 ist von 
25 dem Volumen abhangig, in dem die Phasenumwandlungen erfolgen. 
Es gilt also 

5£l=IiAl5x (11) 

30 P ist dabei eine vorab ermittelbare Proportionalitatskon- 
stante. 

Es ergibt sich ferner eine Anderung der Flache einer Grenz- 
schicht 16 zwischen Ferrit und Zementit. Auch die Anderung 
35 dieser Flache liefert einen Beitrag SE2 zur Energiebilanz. 
Dieser Beitrag 5e2 ist proportional zur Anderung der Flache 
der Grenzschicht 16. Es gilt also 
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6£2 = 2'jA6jc (12) 

7 ist dabei wieder eine vorab erinittelbare Proportionalitats- 
konstante, Der Faktor 2 ergibt sich dadurch, dass pro Phasen- 
grenze 11 r 12 bzw. pro Schicht 13/ 14 je eine Grenzschicht 16 
vorhanden ist und zwei Phasengrenzen 11 , 12 bzw. zwei Schich- 
ten 13, 14 betrachtet warden, 

Anhand der Gleichiangen 11 und 12 lasst sich somit ein kriti- 
scher Lamellenabstand 1' ermitteln, bei dem die Energiebi- 
lanz, welche die beiden Beitrage 5e1 und 8e2 berucksichtigt, 
ausgeglichen ist- Fiir diesen kritischen Lamellenabstand 1' 
gilt 

r=2y/p. (13) 

Der Lamellenabstand 1, also die Summe der Schicht dicken 11 
und 12, wird nun auf das ca. eineinhalb- bis dreifache dieses 
kritischen Lamellenabstand 1' gesetzt, z. B. auf das ca* 
zweifache. Die Schichtdicken 11 bzw. 12 ergeben sich dann zu 

/l=a-?)/ und (14) 

I2 = ql (15) 

Die Schichtdicken 11, 12 sind proportional zu Flachenmafien 
Fl, F2, welche die Phasengrenzen 11, 12 aufweisen. Je eines 
der Flachenmafie Fl, F2 ist in FIG 14 durch Schraffieren ent- 
sprechend hervorgehoben . 

Wenn in einem der Schritte S29, S30 eine Verschiebung 5x' , 
hyi" einer Phasengrenze 11, 12, ermittelt wurde, hat sich so- 
mit die Lage der Phasengrenze 11, 12 geandert. In diesem Fall 
wird in einem Schritt S37 (selbstverstandlich unter Beachtung 
des Vorzeichens der mittleren Verschiebung 8x) eine neue Lage 
der Phasengrenzen 11, 12 oder 11 und 12 bestimmt. In einem 
Schritt S38 wird dann anhand einer nicht linearen Funktion 
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der Lage der Phasengrenze 11, 12 bzw. der Phasengrenzen 11, 
12 der Anteil p3 des Austenits ermittelt. Die nicht lineare 
Funktion beriicksichtigt dabei insbesondere den Umstand, dass 
das Stefan-Problem des Schrittes S29 bzw. die Stef an-Probleme 
5 des Schrittes S30 eindimensional angesetzt und gelost wurden, 
in der Realitat aber eine dreidimensionale Umwandlung er- 
folgt. 

In einem Schritt S3 9 werden sodann die Anderungen der Anteile 
10 pi/ p2 der beiden anderen Phasen Ferrit und Zementit be- 

stimint. Falls dabei die Schritte S37 bis S39 vom Schritt S29 
aus erreicht werden, ist selbstverstandlich der Anteil q Null 
Oder Eins. 

15 Mit der mittleren Verschiebung 5x der Phasengrenzen 11, 12 

und dem Anteil q des Zementit an der Schichtstruktur (0 < q < 
1) lasst sich also auf einfache Weise ermitteln, welche Ande- 
rungen sich fiir die 2Uiteile pi, p2, p3 der Phasen des Stahls 
beziiglich des betrachteten Vol\amenelements 9 ergeben. 

20 

• Trotz der vereinf achenden Annahmen des obigen Modellierver- 
fahrens ergibt sich ein erheblicher Rechenaufwand* Daher wer- 
den gemaB FIG 8 vor dem Schritt S22 Schritte S40 und S41 und 
nach dem Schritt S39 ein Schritt S42 eingefiigt- 

25 

Im Schritt 340 werden die Volumenelemente 9 zu Gruppen zusam- 
men gefasst. Beispielsweise kann eine Zusammenfassung mehre- 
rer Volumenelemente 9 erfolgen, die in Richtung der Band- 
. breite b, der Banddicke d und/oder der Bandgeschwindigkeit v 

30 aneinander angrenzen, Auch Kombinationen hiervon sind mog- 
lich, Im Schritt S41 wird dann pro Gruppe von Volumenelemen- 
ten 9 je ein einziges der Volumenelemente 9 ausgewahlt. Nur 
fur die ausgewahlten der Volumenelemente 9 werden in den 
Schritten S22 bis S39 die Dif ferentialgleichung, das Stefan- 

35 Problem bzw. die Stef an-Probleme gelost und die Anteile pi 
p2, p3 der Phasen rechnerisch ermittelt. 
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GemaB FIG 8 wird bei jeder Gruppe das L6s\Higsergebnis des be- 
trachteten Volumenelements 9 im Schritt S42 fur die anderen 
Volumenelemente 9 der je-weiligen Gruppe ubernoimnen . Dies 
stellt die einfachste Vorgehensweise dar. Es wSre aber auch 
5 denkbar, die Anteile pi, p2, p3 der Phasen von Volumenelemen- 
ten 9, deren Phasenverteilung nicht explizit berechnet worden 
ist, durch lineare Oder nicht lineare Interpolation zu ermit- 
teln. 

10 Auf Grund dieses Umstands, namlich dass bei jeder Gruppe von 
Volximenelementen 9 nur die Konzentrationsverteiiung K eines 
der Volumenelemente 9 erreclinet wird und fiir das L5sen der 
Warmeleitungsgleichung nicht die Konzentrationsverteilung K, 
sondern nur die Anteile pi, p2, p3 der vorhandenen Phasen be- 

15 notigt werden, kann gegebenenf alls auch die Konzentrations- 
verteilung K nur fiir dieses eine Volumenelement 9 jeder 
Gruppe vorgegeben werden. 

Die warmeleitungsgleichung hingegen wird in einem Schritt S43 
20 individuell fiir jedes der Volumenelemente 9 gelost. Der 

Schritt 843 kann dabei alternativ vor Oder nach dem Losen der 
Phasenumwandlungsgleichung in den Schritten S22 bis S42 aus- 
gefuhrt werden. 

25 Das Losen der Warmeleitungsgleichung ist ohne Weiteres mog- 

lich. Denn die Temperaturen T der einzelnen Volumenelemente 9 
sind - siehe obige Erlauterungen zum Schritt S22 - ohne Wei- 
teres ermittelbar, so dass auch deren Ortsgradient sofort er- 
mittelbar ist. Da auch die Anteile pi, p2, p3 der Phasen der 

30 einzelnen Voliimenelemente 9 ebenfalls bekannt sind, lasst 

sich auch die Warmeleitfahigkeit A, der einzelnen Volumenele- 
mente 9 ohne Weiteres ermitteln. Da schliefilich die Dichte p 
im Wesentlichen eine Konstante ist und die Enthalpie H der 
einzelnen Voliimenelemente 9 direkt gegeben ist, ist somit 

35 auch die Warmeleitungsgleichung als Ganzes losbar. 
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Mittels der vorliegenden Erfindung wird also ein physikali- 
schen Modell erstellt, bei dem auf Basis der Gibbsschen 
freien Enthalpie G und des Dif fusionsgesetzes (Stefan-Prob- 
lem) das Umwandlungsgeschehen hinsichtlich entstehender Pha- 
5 senanteile pi, p2, p3 und Umwandliongsgeschwindigkeit mit sehr 
hoher Genauigkeit beschreibbar ist- Es werden allgemein giil- 
tige Ergebnisse ermoglicht, welche auch die Behandlung noch 
unbekannter Materialien und Werkstoffe erlauben. Die Ansatze 
sind dariiber hinaus nicht nur zur Temperaturberechnung/ son- 
10 dern auch zur Berechnung von Gefiigestruktur und KorngroBe an- 
wendbar , 

Das erfindungsgemaBe Modellierverf ahren kann an beliebiger 
Stelle eingesetzt werden, also insbesondere beispielsweise 
15 auch bei Kiihlvorgangen zwischen Walzgeriisten einer WalzstraBe 
Oder bei sogenanntem ferritischem Walzen. Es eignet sich auch 
zur Beschreibung der Riickumwandlung beim Aufheizen des Stah- 
les . 
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Patentanspruche 

1. Rechnergestiitztes Modellierverf ahren fur das Verhalten ei- 
nes Stahlvoliunens (1) mit einer Volumenoberflache, 
5 - wobei ein Rechner (4) anhand eines momentanen Anfangszu- 

stands (ZA) des Stahlvolumens (1) und mindestens einer iober 
die Volxomenoberflache auf das Stahlvolumen (1) einwirkenden 
momentanen Einf lussgrSBe (W) durch L5sen einer WSrmelei- 
tungsgleichung und einer Phasenumwandlungsgleichung einen 
10 Folgezustand (ZF) des Stahlvolumens (1) ermittelt, 

- wobei die mindestens eine Einf lussgroBe (W) fur eine Anzahl 
von Flachenelementen (10) der Volumenoberf ISche jeweils ei- 
nen lokalen Einfluss umfasst und die lokalen Einfliisse uber 
das jeweilige Flachenelement (10) auf das Stahlvolumen (1) 

15 einwirken, 

- wobei der Anf angszustand (ZA) und der Folgezustand (ZF) fiir 
eine Anzahl von Voliomenelementen (9) des Stahlvolxamens (1) 
jeweils lokale Anteile (pl^ p2, p3) von modellierten Phasen 
des Stahls und eine einen lokalen Energieinhalt des Stahls 

20 beschreibende GroBe (H) umfassen, 

- wobei die modellierten Phasen des Stahls Austenit und eine 
erste weitere Phase umfassen, in die Austenit umwandelbar 
ist und die in Austenit umwandelbar ist, 

dadurch gekennzeichnet, 
25 - dass der Anf angszustand (ZA) und der Folgezustand (ZF) fur 
mindestens eines der Volumenelemente (9) auch eine lokale 
Konzentrationsverteilung (K) mindestens eines im Stahl be- 
weglichen Legierungselements umfassen, 

- dass im Rahmen der Umwandlungsgleichung fur das mindestens 
30 eine Volumenelement (9) ermittelt wird, welche Konzentrati- 

onen (kl^ kS; k2, k4) des mindestens einen beweglichen Le- 
gierungselements beidseits einer ersten Phasengrenze (11, 
12) zwischen Austenit und der ersten weiteren Phase vorlie- 
gen, 

35 - dass durch Losen eines ersten Stefan-Problems ermittelt 

wird, ob und wie sich die Konzentrationsverteilung (K) des 
mindestens einen beweglichen Legierungselements im austeni- 
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tischen Bereich des betrachteten Volvunenelements (9) andert 
und ob und um welches AusmaB (8x, 8x' , 6x") sich die erste 
Phasengrenze (11, 12) dadurch verschiebt, und 

- dass die lokalen Anteile (pi, p2, p3) der Phasen anhand ei~ 
5 ner durch das Ausmafi (5x) der Verschiebung der ersten Pha- 
sengrenze (11, 12) bestimmten Lage der ersten Phasengrenze 
(11, 12) ermittelt werden. 

2. Modellierverfahren nach Anspruch 1, 

10 dadurch gekennzeichnet, 

- dass die modellierten Phasen des Stahls auch eine zweite 
weitere Phase umfassen, in die Austenit umwandelbar ist und 
die in Austenit umwandelbar ist, 

- dass fiir das betrachtete Volumenelement (9) im Rahmen der 
15 Dmwandlungsgleichung auch ermittelt wird, welche Konzentra- 

tionen (k2, k4; kl, k3) des itiindestens einen beweglichen 
Legierungselements beidseits einer zweiten Phasengrenze 
(12, 11) zwischen Austenit und der zweiten weiteren Phase 
vorliegen, 

20 - dass durch zusatzliches Losen eines zweiten Stefan -Problems 
ermittelt wird, ob und wie sich die Konzentrationsvertei- 
lung (K) des mindestens einen beweglichen Legierungsele- 
ments im austenitischen Bereich des betrachteten Volumen- 
elements (9) andert und ob und um welches AusmaB (Sx", 6x' ) 

25 sich die zweite Phasengrenze (12, 11) dadurch verschiebt, 

- dass die Stef an-Probleme miteinander gekoppelt sind, 

- dass den Phasengrenzen (11, 12) FlachenmaBe (Fl, F2) zuge- 
ordnet sind, 

- dass ein Anteil (q) des der zweiten Phasengrenze (12) zuge- 
30 ordneten FlachenmaBes (F2) an der Siamme der FlachenmaBe 

(Fl, F2) bestimmt wird und 

- dass die lokalen Anteile (pi, p2, p3) auch vom Anteil (q) 
des der zweiten Phasengrenze (12) zugeordneten Flachenmafies 
(F2) an der Summe der FlachenmaBe (Fl, F2) abhangen, 

35 

3. Modellierverfahren nach Anspruch 2, 

dadurch gekennzeichnet, dass der An- 
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teil (q) des der zweiten Phasengrenze (12) zugeordneten Fla- 
chenmaBes (F2) an der Sxamme der Flachenmafie (Fl, F2) derart 
bestimrat wird, dass die Phasengrenzen (11, 12) stets neben- 
einander angeordnet bleiben, 

5 

4. Modellierverf ahren nach Anspruch 2, 

dadurch gekennzeichnet, dass der An- 
teil (q) des der zweiten Phasengrenze (12) zugeordneten Fla- 
chenmaBes (F2) an der Summe der FlachenmaBe (Fl, F2) derart 
10 nachgefuhrt ward/ dass die Phasengrenzen (11 , 12) aufeinander 
zu streben, 

5. Modellierverf ahren nach Anspruch 2, 3, oder 4, 
dadurch gekennzeichnet, dass anhand 

15 des Anteils (q) des der zweiten Phasengrenze (12) zugeordne- 
. ten FiachenmaBes (F2) an der Suittme der FlachenmaBe (Fl, F2) 
abgeleitet wird/ ob Austenit nur in die erste weitere Phase , 
nur in die zweite weitere Phase oder sowohl in die erste als 
auch in die zweite weitere Phase umgewandelt wird. 

20 

6. Modellierverf ahren nach einem der obigen Anspruche, 
dadurch gekennzeichnet, 

- dass das betrachtete Volumenelement (9) quaderformig ausge- 
bildet ist und drei Grundabmessungen (A, B, C) aufweist, 

25 - dass die erste Phasengrenze (11, 12) als Rechteck mit einer 
ersten Langsseite und einer ersten Querseite ausgebildet 
ist und 

- dass die erste Langsseite mit einer ersten der Grundabmes- 
sungen (A/ B, C) korrespondiert/ die erste Querseite paral- 

30 lei zu einer zweiten der Grundabmessungen (A, B, C) ver- 

lauft und Verschiebungen (5x' , 5x") der ersten Phasengrenze 
(11, 12) parallel zu der dritten der Grundaibmessungen (A, 
B, C) erfolgen. 

'35 7. Modellierverf ahren nach Ansp2ruch 6 und einem der Anspriiche 
2 bis 5, dadurch gekennzeichnet. 
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- dass die zweite Phasengrenze (12, 11) als Rechteck mit ei- 
ner zweiten Langsseite und einer zweiten Querseite ausge- 
bildet ist \and 

- dass die zweite Langsseite mit der ersten der Grundabmes- 
5 sungen (A, B, C) korrespondiertr die zweite Querseite pa- 

rallel zur zweiten der Griindabmes sungen (A, C) verlauft 
und Verschiebungen (Sx", 8xM der zweiten Phasengrenze (12, 
11) parallel zur dritten der Grundabmes sungen (A, B, C) er- 
folgen, 

10 

8. Modellieirverf ahren nach Anspruch 1, 

dadurch gekennzeichnet, dass die 
Summe (1) der Querseiten der Phasengrenzen (11, 12) ungefahr 
gleich dem 1,5- bis 3-fachen eines kritischen Lamellenabstan- 
15 des (1' ) ist, bei dem eine Energiebilanz , die einerseits mit 
dem Verschieben der Phasengrenzen (11, 12) korrespondierende 
Phasenurawandlungen des Stahls und andererseits mit dem Ver- 
schieben der Phasengrenzen (11, 12) korrespondierende Ande- 

■ 

rungen der Flache einer Grenzschicht (16) zwischen der ersten 
20 und der zweiten weiteren Phase beriicksichtigt, ausgeglichen 
ist . 

9. Modellierverfahren nach einem der obigen Anspriiche, 
dadurch gekennzeichnet, dass das Ste- 

25 fan-Problem eindimensional angesetzt und gelost wird bzw. die 
Stef an-Probleme eindimensional angesetzt und gelost werden 
und dass der Anteil (p3) an Austenit anhand einer nicht line- 
aren Funktion der Lage der Phasengrenze (11, 12) bzw. der 
Phasengrenzen (11, 12) ermittelt wird. 

30 

10. Modellierverfahren nach einem der obigen Anspriiche, 
dadurch gekennzeichnet, dass die Er- 
mittlung, in welchen Konzentrationen (kl bis k4) das mindes- 
tens eine bewegliche Legierungselement beidseits der ersten 

35 Phasengrenze (11, 12) bzw, beidseits der ersten und beidseits 
der zweiten Phasengrenze (11, 12) vorliegt, anhand der 
Gibbsschen freien Enthalpien (Gl, G2, G3) der Phasen erfolgt. 
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11. Modellierverf ahren nach einem der obigen Anspruche, 
dadurch gekennzeichnet, dass anhand 
der im Anf angszustand (ZA) bereits vorhandenen Phasen xand an- 
5 hand der Gibbsschen freien Enthalpien (Gl^ G2, G3) der Phasen 
ermittelt wird/ ob sowohl Austenit als auch die erste weitere 
Phase vorliegen bzw. ob zusatzlich zu Austenit und der erst en 
weiteren Phase auch die zweite weitere Phase vorliegt. 

10 12. Modellieirverf ahren nach einem der obigen Anspriiche, 
dadurch gekennzeichnet, dass das 
Stahlvolumen (1) eine Vielzahl von VolTimenelementen (9) um- 
fasst, dass das Stefan -Problem bzw. die Stef an-Probleme nur 
fur einen Teil der Volumenelemente (9) gelost wird bzw. wer- 

15 den und dass die lokalen Anteile (pi, p2, p3) der Phasen der 
anderen Volumenelemente (9) anhand der lokalen Anteile (pi, 
p2, p3) der Phasen des Teils der Volumenelemente (9) ermit- 
telt werden. 

20 13. Modellierverf ahren nach Anspruch 12, 

dadurch gekennzeichnet, dass die 
Warmeleitungsgleichung individuell fur jedes Volumenelement 
(9) gel6st wird. 

25 14. Modellierverf ahren nach einem der Anspriiche 1 bis 13, 
dadurch gekennzeichnet, 

- dass dem Rechner (4) ein Erstzustand (Z) und zumindest eine 
gewiinschte EndgrSBe (f*) vorgegeben werden, 

- dass das Modellierverfahren nach einem der Anspriiche 1 bis 
30 10 iterativ angewendet wird, 

- dass der Anf angszustcind (ZA) der ersten Iteration dem Erst- 
zustand (Z) und der Anfangszustand (ZA) jeder weiteren Ite- 
ration dem unmittelbar zuvor ermittelt en Folgezustand (ZF) 
entspricht und 

35 - dass anhand des nach einer letzten Iteration ermittelten 
Folgezustands (ZF) eine erwartete EndgroBe (f ' ) ermittelt 
und itiit der gewiinschten EndgroJJe (f*) verglichen wird. 
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15. Modellierverfahren nach Anspruch 14, 

dadurch gekennzeichnet, dass es on- 
line und in Echtzeit Oder offline ausgefiihrt wird. 

5 

16. Modellierverfahren nach Anspruch 14 oder 15, 
dadurch gekennzeichnet, dass die Ein- 
flussgroBen (W) der Iterationen in ihrer Gesamtheit einem 
EinflussgroBenverlauf entsprechen, dass der Rechner (4) auf- 

10 grund des Vergleichs der erwarteten EndgroiBe (f) itiit der ge- 
wunschten EndgroBe (f*) den EinflussgroBenverlauf variiert 
und ausgehend vom Erstzustand (Z) das Modellierverfahren nach 
Anspruch 12 erneut ausfuhrt, bis zximindest die erwartete End- 
groBe (f) mit der gewiinschten EndgroJie (f*) korrespondiert . 

15 

17- Modellierverfahren nach einem der Anspriiche 1 bis 13, 
dadurch gekennzeichnet, dass es on- 
line und in Echtzeit ausgefiihrt wird, dass der Rechner (4) 
anhand einer aus dem Anf angszustand (ZA) ermittelten Anfangs- 
20 grofie (f) und einer gewiinschten FolgegroBe (f*) die Einfluss- 
groBe (W) ermittelt und dass der Rechner (4) eine Beeinflus- 
sungseinrichtung (2) derart ansteuert, dass das Stahlvolumen 
(1) entsprechend der ermittelten EinflussgroBe (W) beein- 
flusst wird. 

25 

18. Datentrager mit einem auf dem Datentrager gespeicherten 
Computerprogramm (6) zur Durchfiihrung eines Modellierverf ah- 
rens nach einem der obigen Anspriiche. 

30 19. Rechner mit einem Massenspeicher (8), in dem ein Compu- 
terprogramm (6) hinterlegt ist, so dass bei Aufruf des Compu- 
terprogramms (6) von dem Rechner ein Modellierverfahren nach 
einem der Anspriiche 1 bis 17 ausfiihrbar ist. 

35 20. Beeinflussungseinrichtung zum Beeinflussen der Temperatur 
eines Steihlvolumens (1) , insbesondere Kiihlstrecke, 
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dadurch gekennzeichnet/ dass sie von 
einem Rechner (4) nach Anspruch 19 gesteuert wird. 

21. Stahl, 

dadurch gekennzeichnet^ dass er eine 
Beeinf lussungseinrichtung (2) nach iVnspruch 20 durchlaufen 
hat, wobei die Einflussgrofie (W) nach Anspruch 16 oder 17 er- 
mittelt wurde. 
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FIG 6 
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FIG 18 

\ 

lose K't-DK'xx=L1 
und KVDK"xx=L2 

mit L1 = -2^f^?j^-qt(1-q)(3q+1) 
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bestimme q derail, dass 
gilt: 6x'=5x" 
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